VSX1 gene analysis in keratoconus.

PURPOSE
To screen the visual system homeobox 1 (VSX1) gene in keratoconus patients.


METHODS
The enntire coding region of VSX1, including intron-exon boundaries were amplified in keratoconus cases (n=50) and controls (n=50). All sequences were analyzed against the ensemble sequence (ENSG00000100987) for VXS1.


RESULTS
Sequencing analysis showed four alterations (p.A182A, p.R217H, p.P237P, and g.25059612C>T) in VSX1 of which g.25059612C>T (in intron 2) was found to be novel. Of these four, p.A182A and p.P237P were present in both cases as well as controls while p.R217H and g.25059612C>T were limited to cases only. All these changes were non-pathogenic.


CONCLUSIONS
In our study no pathogenic VSX1 mutation was identified. The role of VSX1 in the pathogenesis of keratoconus is still controversial. VSX1 mutations are responsible for a very small fraction of all observed keratoconus cases. The absence of pathogenic mutations in VSX1 in our patients indicates that other genetic loci like 13q32 as suggested by a recent study may be involved in the pathogenesis of this disorder.

for screening mutations in posterior polymorphous corneal dystrophy (PPCD) and keratoconus [8]. VSX1 is considered important in ocular development and is particularly involved in the developing cornea. Expression in human was demonstrated in embryonic craniofacial, adult corneal, and adult retinal cDNA libraries [16]. VSX1 mRNA has been found in the outer tier of the inner nuclear layer of the human retina and the cornea [17]. Mutations in this gene are also associated with posterior polymorphous dystrophy. VSX1 is highly conserved across many species [18][19][20]. Several mutations, such as p.D144E, p.G160D, p.P247R, p.L159M, p.R166W, and p.H244R have been reported by various groups [14,17,21] but a definite pathogenic role of these mutations in keratoconus is not yet established. In this study we present the results of VSX1 gene analysis in 50 keratoconus patients and controls from north India.  videokeratography (VKG), Orbscan, visual testing, fundoscopy, slitlamp-biomicroscopy, and retinoscopy. Of these patients, 29 were males and 21 were females. The mean age of presentation was 18.2 years. Diagnosis of keratoconus involved the presence of characteristic topographic features, such as inferior or central corneal steepening, or an asymmetric bowtie pattern with skewing of the radial axes, and the presence of one or more of the following characteristic clinical features in one or both eyes: conical corneal deformation, munsen sign, corneal stromal thinning, a Fleischer ring or Vogt striae. All cases were sporadic without any family history. All keratoconus cases secondary to causes like trauma, surgery, Ehlers Danlos syndrome, Osteogenesis Imperfecta, and pellucid marginal degeneration were excluded from the study.
After informed consent, detailed personal, medical and occupational history was collected and a family tree up to three generations was drawn. Fifty ethnically matched normal individuals without any ocular disorder were enrolled as controls. Health information was obtained from controls through the questionnaire; all underwent ophthalmological examination and blood sample (5 ml) was collected in EDTA (EDTA) vacutainers (Greiner Bio-One GmbH, Frickenhausen, Germany) from patients and controls for DNA extraction. DNA was extracted from whole blood samples of all patients and controls using the phenol-chloroform method. PCR and DNA sequencing: All the coding regions of VSX1 including exon-intron junctions were amplified using a set of five oligonucleotide primers (Table 2) All PCR products were analyzed on 1.8% agarose gel stained with ethidium-bromide (EtBr; 10 mg/ml). Agarose gels were analyzed using gel documentation system (Applied Biosystems, Carlsbad, CA). Amplified PCR products were purified using gel/PCR DNA fragments extraction kit (DF100; Geneaid Biotech Ltd., Sijhih City Taiwan). The purified PCR products were sent for sequencing at MCLAB (Molecular Cloning Laboratories, South San Francisco, CA) Nucleotide sequences were compared with the VSX1 ensembl reference sequence.
Insilico analysis of missense mutations: Two homology based programs PolyPhen (Polymorphism Phenotyping) and SIFT (Sorting Intolerant From Tolerant) analysis tool were used to predict the functional impact of missense changes identified in this study.
PolyPhen structurally analyzes an amino acid polymorphism and predicts whether that amino acid change is likely to be deleterious to protein function [22,23]. The prediction is based on the position-specific independent counts (PSIC) score derived from multiple sequence alignments of observations. PolyPhen scores of >2.0 indicate the polymorphism is probably damaging to protein function. Scores of 1.5-2.0 are possibly damaging, and scores of <1.5 are likely benign.
SIFT is a sequence homology-based tool that sorts intolerant from tolerant amino acid substitutions and predicts whether an amino acid substitution in a protein will have a phenotypic effect [24,25]. SIFT is based on the premise that protein evolution is correlated with protein function. Positions important for function should be conserved in an alignment of the protein family, whereas unimportant positions should appear diverse in an alignment. Positions with normalized probabilities less than 0.05 are predicted to be deleterious and, those greater than or equal to 0.05 are predicted to be tolerated. We have also used improved Splice Site predictor tool [26] for prediction the effect of an intronic nucleotide change on splicing.

RESULTS
DNA sequencing analysis of 50 patients and 50 controls revealed a total of four nucleotide changes (Table 3)   one was novel and 3 have been previously reported. Details of these changes are given below.
Arginine217Histidine (p.R217H): In this mutation a single nucleotide thymine (T) was replaced by adenine at  position g.25059442 (rs6138482); cDNA position c.650; codon 217. This change resulted in a codon change from CGC>CAC resulting in non-synonymous change p.arg217his (p.R217H; Figure 2) in protein. This change was present in only one case and was homozygous but was absent in controls.
Proline237Proline (p.P237P): In this mutation a single nucleotide T was replaced by A at position g.25059381 (rs56157240); c.711; codon 237 resulted in a codon change CCT>CCA which predicts a synonymous change p.pro237pro (p.P237P; Figure 3). This change was present in 18/50 cases (7 were homozygous and 11 were heterozygous) being also present in 14/50 controls (9 were homozygous and 5 were heterozygous).
Cytosine to Thymine in intron 2: A novel single nucleotide change C>T at g.25059612 (Figure 4) was present in three cases but absent in controls. Alteration is located in 2nd intron (IVS3-24C>T). This change was registered at GenBank with accession number GU471016.
PolyPhen and SIFT analysis of p.R217H showed that it is non-pathogenic (SIFT score is >0.05 and PSIC score is <1.5).
Improved splice site predictor tool analysis of g. 25059612C>T showed that this location (g.25059612) is not present at splice site and may not create splicing error in VSX1 mRNA. Conclusively, no pathogenic change was identified in our patients.

DISCUSSION
In this study we analyzed VSX1 in 50 unrelated keratoconus patients and controls from north India. In our patient's, males were affected more than females. Mutations in VSX1 gene have been identified in association with keratoconus [12,14,15,17,21]. Human VSX1 is a member of the CVC domain containing paired-like class of homeoproteins. VSX1 expression in humans was detected in embryonic craniofacial, adult retinal, and adult corneal tissues [17,27]. The role of VSX1 in keratoconus is still ambiguous. Previous studies have shown that the pathogenesis of keratoconus is very complex and several genes and gene environmental interactions play a critical role in disease prognosis. In fact, VSX1 may have a pleiotropic action among the tissues of the cornea leading to keratoconus in some cases as observed for the transforming growth factor (TGFBI) gene; which causes four distinct autosomal dominant corneal diseases [28].
In this study, four sequence variations were detected of which 3 have been previously reported [15,20] and one was novel. SIFT and PolyPhen analysis of p.217H showed that it to be a non-pathogenic change. Similar findings have already been reported in European population [15].
The VSXI variants reported in various studies include p.L17P, p.D144E, p.N151S, p.L159M, p.G160V, p.G160D, p.R166W, p.Q175H, p.H244R, and p.P247R [12][13][14][15]17,21]. Some of these were initially reported to be pathogenic but their pathogenicity could not be confirmed as some of these variants were observed in unaffected individuals also. He'on and associates [14] identified a compound heterozygous change with p.P247R and p.G160D and reported p.G160D to be pathogenic and p.P247R to be nonpathogenic. Another study reported p.P247R as a pathogenic change because it was found to be co-segregating with keratoconus. Similarly, p.D144E mutation was initially reported as pathogenic [20, 21,] but subsequent studies identified its presence in unaffected individuals and suggested this to be a nonpathogenic polymorphism [29,30]. The variants p.R166W, p.H244R, and p.L159M have been identified in keratoconus patients but these changes did not segregate with the disease phenotype in their family members and hence were not considered sufficiently significant to support a pathogenic role in keratoconus [31]. Similarly, variants p.G160V and p.N151S have been identified in patients from the Korean population [12] but these changes have not been reported in other populations.
In a recent study (2009) from India, VSX1 was screened in 66 keratoconus cases and a potentially pathogenic change (p.Q175H) was identified in one case only [13]. In this study, the second from India, no pathogenic change was indentified in VSX1. Similar results have been published recently [29,32,33]. So lack of possibly pathogenic changes in VSX1 gene in keratoconus patients suggests that mutations of VSX1 could only be responsible for a very small fraction of all observed cases and need to be investigated in different populations. This also suggests that other genetic loci like 13q32 as suggested by Gajecka et al. [34] may be involved in the pathogenesis of keratoconus.